Focusing through random media: eigenchannel participation number and intensity 
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Using random matrix calculations, we show that, the contrast between maximally focused intensity 
through random media and the background of the transmitted speckle pattern for diffusive waves is, 
Miv = 1 + Neff, where N^ff is the eigenchannel participation number for the transmission matrix. 
For diffusive waves, N^ff is the inverse of the degree of intensity correlation, k. The profile of 
the focused beam relative to the ensemble average intensity is expressed in terms of the square of 
the normalized spatial field correlation function, F{Ar), and k. These results are demonstrated in 
microwaves experiments and provide the parameters for optimal focusing and the limits of imaging. 

PACS numbers: 42.25.Dd, 42.25.Bs 



Electromagnetic, acoustic and electronic waves are 
used to transmit energy and information and to control, 
probe and image our environment. The particularities of 
wave interactions in diverse systems is the basis of our 
rich experience of the world, but their common charac- 
teristics provide a framework for understanding transport 
and often point to new applications. 

Waves traversing a random medium are typically 
so completely scrambled that they produce a random 
speckle pattern of transmitted flux bearing no relation 
to the incident waveform. Nonetheless, intensity can be 
focused at any point by phasing all sources of the wave 
so that the fields arriving at the selected point from these 
sources add constructively. Light has recently been fo- 
cused through a random medium by employing a genetic 
algorithm to manipulate the incident field in each gen- 
eration based on feedback of transmitted intensity at a 
pointi, and by recording the transmission matrix and us- 
ing phase conjugation.^ Light can also be focused through 
a random slab in both space and time by manipulating 
the incident fieldi^""— 

Maximum focused intensity is achieved by construct- 
ing an incident field by phase conjugating the Green's 
function from the focal point to the incident surfacci^ 
Phase conjugation of monochromatic waves is analogous 
to time reversal for pulsed signals which is widely used 
to focus ultrasound and electromagnetic radiation within 
reverberant cavitiesj^i^ or through random media^*^'^^. 

The contrast and spatial profile of focused intensity in 
transmission through a multiple scattering sample must 
depend upon the random medium as well as upon the 
incident field. It is natural to explore the role of the 
medium via the complete set of the A'^ orthogonal eigen- 
channels of the transmission matrix and their associated 
transmission eigenvalues t„. The field transmission ma- 
trix t relates the fields Ea and Ef, between incoming chan- 
nel a and outgoing channel b with, Ei, — "YlH ^baEa- 

The full transmission matrix is not generally available; 
however, focusing through random media may also be 
characterized via the correlation and the fiuctuations of 
the transmitted intensity, which are accessible experi- 
mentally. These statistical characteristics of the inter- 



action of the wave with the medium reflect the number 
of statistically independent constituents of the transmit- 
ted field. Intensity correlation across the transmitted 
speckle pattern persists beyond the range over which the 
field is correlatedi2r2& and may be characterized through 
the degree of intensity correlation This is the value 
of the correlation of intensity normalized by its ensem- 
ble average between points at which field correlation 
vanishes. Correlation in intensity and fluctuations of 
total transmission normalized by its ensemble average, 
Sa = Ta/{Ta), whcrc Ta = Yl,b l^f)a|^, are linked via the 
equality, k — var(sa) for N ^ 1.^^ k is also inversely pro- 
portional to the dimensionless conductance, g ji3ii4,i6- _20 
which is the ensemble average of the optical transmit- 
tance, g = (T) = (E^, Ihaf) = rn)-'' 

Though focusing through random media has been 
demonstrated, the focused intensity pattern has not been 
fully characterized and related to experimentally accessi- 
ble parameters. In addition, the full transmission matrix, 
which reflects the correlation within the medium has not 
been measured and related to focusing. In the present 
paper we show that, for maximum focused intensity, the 
ensemble average of the contrast in focusing may be ex- 
pressed in terms of the degree of intensity correlation, k, 
and the number of effective eigenvalues, N^ff, which is 
the eigenchannel participation number of the transmis- 
sion matrix. We show that in the diffusive limit, the 
average of intensity at a distance Ar from the focus rela- 
tive to the ensemble average of the intensity is given by, 
{Ifoc{Ar))/{I} = N [(FiAr) + k) / (1 + k)]. Here i^(Ar) 
is the square of the field correlation function with dis- 
placement normalized to unity at Ar = and k is equal 
to the inverse of N^f / . In the limit of strong localization, 
the contrast approaches unity and the wave can no longer 
be focused. These results are confirmed via microwave 
measurements of the focusing profile and contrast, Neff, 
and K in random ensembles of samples over a wide range 
of Neff. The description of focusing in space in terms 
of the degree of intensity correlation can be extended to 
describe focusing in both time and space. 

Focusing with maximum intensity on a point /3 for a 



2 



normalized incident field is achieved by phase conjugat- 
ing the field transmission coefficient between the target 
point at /3 and input points a so that the incident field is, 

Tp.^ Components of the 



Ea 

transmitted field originating from different points then 
interfere constructively at /3 to give the focused intensity, 
which is equal to the total transmission through the op- 

2 ^ 

posite surface for a source at (3, Ip — X^a l^^^l /Tp = 
Tp. For M incoming channels with M < N, where N 
is the number of independent channels supported in the 
space surrounding the medium, (Ip) is enhanced by a 
factor M over (/), (Ip) = M (I) = Mg/N"^. On the 

other hand, the background intensity at points 6 7^ /3 is, 
2 

/Tp. 

Decomposing the field into orthonormal incoming and 
outgoing singular vectors v„ and u„, respectively, t can 
be written as t = UAV = ^n=i^n^n^li- The singu- 
lar values A„ are the square root of the eigenvalues r„ 
obtained by diagonalizing the Hermitian matrix tt^ and 
the singular vectors v„ and u„ are the waveforms at in- 
put and ouput surfaces, respectively, which couple selec- 
tively to the nth eigenchannel propagating through the 
medium. The focused and background intensities can be 
expressed in terms of the singular vectors (see Appendix 
A), 



= Tn\Unp\'^ 
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Since the background intensity is the result of interfer- 
ence between randomly phased statistically independent 
elements associated with different eigenchannels, the con- 
trast — {Ip/h) depends on the effective number of 
eigenchannels contributing to the transmission. 

The average of the contrast between the peak and 
background intensity for M input points is shown in Ap- 
pendix A to be. 



Ma/ = 



M 



(2) 



where Ne/f = /Er^=i^n), the eigenchan- 

nel participation number of the transmission matrix, is 
the effective number of uncorrelated speckle patterns con- 
tributing to the transmitted field. 

N^ff reflects the internal statistics of t„ and the 
statistics of relative intensity and is independent of ti. 
Dorokhov^^ and Imry^"^ have shown that for good con- 
ductors, g is equal to the number of open or active eigen- 
channels for which r„ > 1/e. This number is propor- 
tional to Neff in nondissipative, diffusive samples. For 
Neff < M, Eq. (2) yields, 



(3) 



Equation (2) is consistent with the previous finding 
that the contrast increases linearly with AI when M ^ 
g.— In this case, the components of an NxAI matrix for 
M incoming channels are uncorrelated and the normal- 
ized singular values of this matrix, A — ^/ \J'^n=i '^n/M 
follow the quarter-circle law for M ^ l.?4 This distri- 
bution gives Neff — M/2 which, taken together with 
Eq. (2), leads to hm — M. The quarter-circle law was 
found in optical^ and acoustio^^ measurements in which 
the number of measured channels M was smaller than g. 

In general, the contrast in focusing depends both upon 
t and on the character of the incident field. The contrast 
increases with the effective number of eigenchannels con- 
tributing to the intensity at the focus. For an arbitrary 
ensemble of incident fields with complex amplitude a„ 
for the ri"^ eigenchannel, the incident field which adds 
coherently at the focal point P is Ea = '^nU^pVna 
and the eigenchannel participation number in the focused 

field becomes N^f^ = ((X^ii |a„| A„) knl^A^). 
When the incident field is constructed by means of phase 
conjugation, a„ — A„ and Nj;^j — N^ff- For the choice, 

a„ = 1/A„, known as the inverse filter^ N^jj = N and 
fiN ^ N. However, since smaller values of A„ are em- 
phasized in this approach, the peak intensity is small. 
The noise level can be high in this case since significant 
weight is given to small values of A„ in which the relative 
contribution of noise is larger. For any set of a„, the con- 
trast achieved when these waves are summed to interfere 
constructively at a point is given by Eq. (2) with Nj^^J 
substituted for Neff. 

Microwaves measurements of the field transmission 
matrix allow us to test the calculations of contrast above 
as well as the relationship between Neff and intensity 
correlation and the spatial distribution for the focused 
wave over a large range of Neff. The comparison is 
made with high resolution for ensembles of statistically 
equivalent disordered samples. Measurements are made 
in samples of randomly positioned alumina spheres con- 
tained in the copper tube in two frequency ranges 31 The 
wave is localized in the frequency range between 10 and 
10.24 GHz and diffusive between 14.7 and 14.94 GHz 
with ^ 30 and 66, respectively. Measurements are 
carried out for sample of lengths i=23, 40 and 61 cm for 
which g ranges from 6.9 to 0.17 and Neff from 13.8 to 
1.16. Measurements are made for horizontal and vertical 
linearly polarized components of the field by rotating the 
source and detector antennas. The antennas are trans- 
lated over a 9 mm grid on the front and back surfaces 
of the waveguide. The field correlation falls rapidly and 
F{Ar = 9mm) 0. The transmission matrix t is com- 
puted using A/2 points for each polarization. New con- 
figurations are created by briefly rotating and vibrating 
the sample tube. 

Typical speckle patterns for diffusive and localized 
waves are displayed in Fig. 1(a) and (b), respectively. 
Patterns produced by focusing the field at the input via 
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FIG. 1. (Color online) Intensity speckle pattern generated for 
L=23 cm for diffusive waves (a) and for L—61 cm for localized 
waves (b) normalized to the average intensity in the respective 
patterns. Focusing at the central point the same frequency as 
in (a) and (b) via phase conjugation is displayed in (c) and 
(d) with 66 and 30 input points, respectively. 



the phase conjugation protocol at a point in the center of 
the output speckle pattern, /3 = (x, y) = (0, 0), are shown 
on Fig. 1(c) and (d). The Whittaker-Shannon sampling 
is used to obtain high resolution patterns.-2^ The focused 
patterns are obtained utilizing the reciprocity of propa- 
gation to calculate the incident field t*^^. Only for dif- 
fusive waves does the focal spot emerge from the back- 
ground. Several closely spaced eigenvalues contribute to 
the transmittance T as can be seen in Fig. 2(a). For 
the diffusive sample of length L=23 cm and g = 6.9, 
we find N^jf = 13.8. In contrast, the first eigenchanncl 
dominates transmission for localized waves so that the in- 
tensity pattern is not very different from the pattern for 
this eigenchanncl. In the sample of length L=61 cm and 
g=0.17, (ti) / (t2) 24 and Nejf = 1.16. As a result, 
it is not possible to significantly enhance the intensity 
at a particular point for localized waves. The intensity 
at the focus is still N (/) but this essentially the average 
intensity in the pattern when only the first transmission 
eigenchanncl is excited. 

Focusing diffusive waves via phase conjugation may 
greatly enhance the background intensity by giving more 
weight to more highly transmitting eigenchannels.- From 
Eq. (1), the total transmission obtained by 

summing the intensity over the output when the wave 
is focused for maximum intensity is given by C4 2 — 
(E„T2) /(E„r„) = g/iVe//,6 instead of (r„) - g/7V 
for illumination by a random wave. This enhances the 
background in the focused pattern over the average in- 
tensity (/) by a factor N/Neff. For the diffusive samples 
studied here with L=23 cm, the total transmission is en- 
hanced five-fold. The values of (r„) are seen in Fig. 2(a) 
to fall exponentially for these samples. This exponential 
falloff corresponds to the distribution, P (t) = g/ r for the 
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FIG. 2. (Color online) (a,b) Averages of r„ for the first 30 
transmission eigenvalues (a) for diffusive waves, L=23 cm, and 
(b) for localized waves, L=61 cm. 



transmission eigenvalues, and give an average value of to- 
tal transmission of g/A^e// = 1/2. Wefindg/A^e// = 0.50 
and 0.53 for L=23 and 61 cm, respectively, giving en- 
hancements of total transmission by factors of 4.8 and 
10 for the focused wave. This results differ from the 
result predicted for the bimodal distribution, in the dif- 
fusive hmit^i^Si^ of g/Neff = 2/3. The departure we 
find from the bimodal distribution may be due to the 
measurements protocol in which the field is measured 
between points on a grid instead of between waveguide 
modes. ''^ In addition, the fractional field correlation has 
a residual value of ~ 0.05 as opposed to the vanishing 
of field correlation between orthogonal channels which 
would occur with continuous spatial wavefronts. 

The contrast hm for M incoming channels is seen in 
Fig. 3 to increase linearly before beginning to saturate. 
The departure from linearity was noted by Derode et al. 
in acoustical studies. Assuming that the average of the 
background intensities is linear in M for M < N gives 
{h) = (/) [1 + {l/^iN ~ l/N) (M - 1)], which leads for 
A > 1 to, 



M 



[l + {l/^iN-l|N) {M-i)Y 



(4) 



Good agreement with measurements is seen in Fig. 3. 

The eigenchannel participation number, N^ff, may be 
directly linked to fluctuations and correlations of inten- 
sity for diffusive waves for quasi-one dimensional sam- 
ples. In such samples, the length greatly exceeds the 
transverse dimensions of the sample with reflecting sides 
so that waves from all incident points are perfectly mixed 
on the sample output. Intensity correlations remain even 
when the field correlation function vanishes. The cumu- 
lant correlation function of normalized intensity versus 
displacement, C (Ar) = {SI{r)6I{r + Ar)), may be ex- 
pressed in terms of the square of the normalized field cor- 
relation function, Fe (Ar) = {E{r)E*{r + Ar)) and n as 
C(Ar) = i^(Ar)-FK(l + F(Ar))J^iilHere, /(r) = \E\^ = 
\tab\\ i{r) = I{r)l{I{r)), E{r) = £;(r)/^vW)) and 
F — F^. K involves both long-range and infinite-range 
correlations. Correlations in intensity and fiuctuations 
of total transmission normalized by its ensemble aver- 
age, Sa = Ta/{Ta), arc linked with k = var (sa), which 
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FIG. 3. (Color online) The contrast is plotted with re- 
spect to the number of input points M for Z/=61 cm (ma- 
genta filled circles), 40 cm (cyan asteriks) and 23 cm (blue 
circles) for difi'usive waves and L=61 cm (red stars) for local- 
ized waves. The curves are the fit of Eq. (4) to the data. 
The black dashed line corresponds to ~ M . In the insert, 
measurements of are compared to (1 + 1/k) for diffusive 
waves with k, equal to the value of k = var(sa) obtained ex- 
perimentally. 
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FIG. 4. (Color online) The ensemble average of normalized 
intensity for focused radiation (blue circles) is compared to 
Eq. (5) (blue solid line) for L=61 cm for diffusive waves. For 
localized waves for L=61 cm, k is replaced by l/(/ijv — 1) in Eq. 
(5). F (Ar) (blue dots) is fit with the theoretical expression 
obtained from the Fourier transform of the specific intensity 
(dashed blue line). The field has been recorded along a line 
with a spacing of 2 mm for 49 input points for L—61 cm. The 
black dashed line is proportional to (J) / {Ifoc (0)) = 1/N. 



is shown in Appendix B to be for diffusive waves and 
iV > 1, 

var (Sa) = K = '^/Neff. (5) 

For the ensemble of diffusive samples of length L—61 
cm, var (sa) = 0.147 is in good agreement with l/N^ff ~ 
0.145. Equations (3) and (5) then give, = 1 + 
1/ var (sa) = 1 + 1/k which is in a good agreement with 
measurements in Fig. 3. 

Significantly higher spatial resolution is achieved by 
focusing in disordered systems than in free space since 
the resolution is not limited by the aperture of the emit- 
ting array. Instead, the resolution is equal to the field 
correlation length which is the inverse of the width of 
the fc- vector distribution of the scattered waves Ji^^Ti^ 
The spatial variation of intensity in the focused speckle 
pattern {Ifoc (A'')) reflects the decay of the ensemble av- 
erage of the coherent sum of eigenchannels at the focus 
towards the average value of the incoherent sum. The 
average intensity pattern for focused waves normalized 
to the value at the focus is derived in Appendix C as, 

{Ifoo (Ar)) ^ < Ell rl> + < Ell T-„>^f (Ar) 
<Elir,^> + <Eli-«>^ 



The ensemble average of the intensity in the focused 
patterns is shown in Fig. 4 and seen to be in an ex- 
cellent agreement with Eq. (6) using measurements of 



the square of the field correlation function F(Ar) and 
At. K, is obtained from the fit of the cumulant intensity 
correlation function C(Ar). i^(Ar) shown in Fig. 4 is 
well fitted with theoretical expressiona^^^ based on the 
equality of i^E(Ar) with the Fourier transform of the 
specific intensity. For localized waves, Eq. (6) no longer 
holds, but good agreement with the focused pattern is 
obtained when k in Eq. (6) is replaced with — 1) 

with obtained experimentally. 

We have considered theoretically focusing in nondis- 
sipative, quasi-one-dimensional samples in steady state. 
The impact of absorption and the nature of focusing in 
the slab geometry, in which the wave is not restricted 
in transverse directions, and in the time domain can be 
outlined in the context of these results. Focusing is lit- 
tle affected by absorption since k is much less sensitive 
to absorption than is transmission.— N^/f is similarly 
insensitive to absorption since the ratio of transmission 
eigenvalues is only slight reduced by absorption. As a 
result of the drop in k and the increase in N^f / , contrast 
in focusing increases slightly with absorption. In a slab 
geometry, k varies with separation between points since 
the overlap of the portion of the medium explored by the 
waves arriving at different points on the output falls with 
separationjiS. 

Equation (6) can be generalized to give the focused 
pattern in space and time. The intensity contrast within 
the focused pattern at a time At relative to the time 
t at which the wave is focused fj,N{t,At) depends on 
both t and At. The intensity in the pattern due to a 
pulse focused at an arbitrary point within the speckle 
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pattern and at a time t is given by replacing F in Eq. 
(6) by F{Ar,At) = F{Ar)F{At),^'^- and k by the 
time varying degree of correlation Ka-{t) for an incident 
pulse with bandwidth a}"^ Here F{At) F|(Ai) is the 
square of the normalized field correlation function with 
time, FE^At). ^^{t) falls with increasing bandwidth as 
a greater number of quasi-normal modes and similarly 
a greater number of transmission eigenchannels are sub- 
stantially excited by the pulse. This gives, 



{Ifoc{Ar,t,At)) F{Ar,At) + K^{t) 



NiHt)) 



(7) 



EM * 



and, 



EM 2 
^ Tn\Unp\ ■ 
n—1 



(Al) 



Whereas the intensity at the focus is the square of a co- 
herent sum, the background intensity reflects the random 
phasing between eigenchannels. For separation much 
greater than the field correlation length, (5r, the aver- 
age background is the sum of the square amplitudes of 
the transmission eigenchannels. 



Ka (t) typically falls for short delay times before increases 
at later times. This has been explained in terms of the 
temporal variation of the number of quasi-normal modes 
of the medium that contribute substantially to transmis- 
sion and, again, the effective number of distinct trans- 
mission eigenchannels. At short times, transmission is 
dominated by a few short-lived modes which surrender 
their energy rapidly, while at later times transmission is 
due to a small number of long-lived modes. As a result 
correlation is high for short and long times but lower at 
intermediate times 4i Finding the transmission matrix at 
different delay times would yield the effective number of 
eigenchannels at a delay time t, Neff(t), and the contrast 
in the speckle pattern accordingly. 

In conclusion, we have provided a description of focus- 
ing through scattering media in space and time for dif- 
fusive and localized waves in terms of the effective eigen- 
channel number and the field correlation function. For 
diffusive waves, the contrast may be expressed in terms of 
an equivalent statistical localization parameter, the de- 
gree of correlation k, which is the inverse eigenchannel 



n: 



In the localized limit, 



participation number, k — ^^eff- 
the contrast approaches unity since transmission is car- 
ried by a single transmission eigenchannel. These results 
are confirmed in microwave measurements and provide a 
systematic framework for focusing radiation for applica- 
tions in imaging and communications. 
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Nos. DMR-0907285 and DMR-0958772 MRI-R2 and by 
the Direction Generale de I'Armement (DGA). 



Appendix A: DERIVATION OF THE CONTRAST 

The contrast fjLM for M incident channels is the 
average of the ratio of the focused intensity, = 



/Tg — Tg, and the background intensity. 



(Ip/Ib). The intensity 

at the focus can be expressed in terms of the trans- 
mission matrix by decomposing the field into singular 
values. The singular vectors are normalized so that 



X„.u 



^M 



/Tp. 



(A2) 



Since the spacing between two points on the measure- 
ment grid is equal to or larger than 6r, we assume that 
components of the singular vector associated with differ- 
ent points on the sample output for the nth eigenchannel, 
Unp and Unb, are uncorrelated. Using Eq. (1), the con- 
trast iiM — {Ip /lb) in the diffusive limit can be expressed 



as. 



.M 



' ^ r. 



(A3) 



For diffusive waves, correlation between the numer- 
ator and the denominator is small. The numerator 
(Ij) can be written as, (/|) = {J2n=iTnWnpf) + 

{J2n=lY.n'^n'TnTri'\Unpf\Un'p\'^). Thc singular VCCtOrS, 

normalized so that (|u„^|^) = 1/M, are independent of 
the eigenvalues. The intensity of the Gaussian singular 
vectors at the output, A/|u„^| , have a negative exponen- 
tial distribution, giving (M^|u„^|'*) — 2. In the diffusive 
limit for systems with time reversal symmetry, var (T) = 
(T2) - {Tf « 2/15,i2ii4i38 while (T) = g is large, so 



that, (r2) « (T)2, and , 



((E^ll^n)")- 



This gives, (/|) « r,?) + rn}') / M\ 

The average denominator in Eq. (A3) is averaged over 



all points b ^ P and is equal to, {J2n=i^nj /-^^ ~ 
{Il)/M\ 

The contrast can finally be expressed as a function 
of the eigenvalues r„, 



(A4) 



Since the eigenchannel participation number. 



{(J2n=i'^n) /E^ii^n). IS equal to 



{J2n=i'^n) /(Z^nli'^n) diffusivc limit,'^ this 

gives Eq. (2). 



^M 



6 



Appendix B: RELATION BETWEEN Neff AND 
var(sa) 



The variance of normalized total transmission, 
Sa = Ta/{Ta), where, = Y.b \*baf ' is 

\2 



var(s,) = {Tl)l{Ta) - 1. We obtain, var(s,) = 



the diflFusive limit, this leads to. 



In 



var (sq) 



(Bl) 



Appendix C: DERIVATION OF THE FOCUSED 
SPATIAL INTENSITY PATTERN 



The numerator on the right hand side of Eq. 
(CI) is also written as, ('^n)('""f'</3<f,u«/3) + 

J2n=i 'En'^n (Tnr„') (u„b<^<,^w„'/5) . Since the compo- 
nents of the singular vectors, Unb, are circular Gaussian 
variables, the average product (unbU^isU'^b'Unis) can be 
broken into the sum, | ^)|^ + {\unbf){\uni3\'^)^^ In 
this expression, |(MnfcU*^)|^ is the square of the field cor- 
relation function of the singular vector Unb, which can be 
approximated as, |(MnfcU^^)|^ — F (Ar) /N'^. This gives, 
{nnbK^u*„,Unp) = F(Ar) /N^ + l/N\ 

For n ^ n' , the singular vectors Unb and Un'b are uncor- 
rected so that (u„6U*^U*,(,U„'/3) = {UnbU*^i3){u*^>bUn'fj) = 

F (Ar) /N"^. The numerator in Eq. (CI) can then be 



written as, ( X/n=i '''nUnbU^p 

) = [(E„=i T-n) F (Ar) 
{En=i'^n)]/N^- Equation (CI) finally leads to. 



N 



When focusing is achieved by phase-conjugation at the 
selected point /3, the average intensity at points b dis- 
placed by Ar from /?, {I foe i^f)) can be expressed as 
a function of the degree of extended correlation k and 
the square of the field correlation function, F (Ar) — 
\{E{r)E*{r + Ar))f/[{I{r)}{I{r + Ar))]. Using Eq. 
(1), the average intensity in the focused pattern normal- 
ized to the average intensity at the focus is. 



N{I) 



N 



^n'^nb^' 



(CI) 



(//cc(Ar)) ^ rn)^F (Ar) + (E 

N{I) 



N 



(C2) 



By virtue of Eq. (Bl) and k — var(sa), Eq. (C2) can 
be rewritten to give the spatial variation of the normal- 
ized intensity pattern with the displacement Ar in terms 
of K, 



{I foe (Ar)) ^ F{Ar) + K 

N{I) 1 + K 



(C3) 
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